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ABSTRACT 

In order to understand the formation and evolution of Magellanic-type dwarf irregular (dim) 
galaxies, one needs to understand their three-dimensional structure. We present measurements 
of the stellar velocity dispersion in NGC 1569, a nearby post-starburst dim galaxy. The stellar 
vertical velocity dispersion, er z , coupled with the maximum rotational velocity derived from H I 
observations, K ma j., gives a measure of how kinematically hot the galaxy is, and, therefore, indi- 
cates its structure. We conclude that the stars in NGC 1569 are in a thick disk with a V maK /a z 
= 2.4 ± 0.7. In addition to the structure, we analyze the ionized gas kinematics from O III 
observations along the morphological major axis. These data show evidence for outflow from 
the inner starburst region and a potential expanding shell near supermassive star cluster (SSC) 
A. When compared to the stellar kinematics, the velocity dispersion of the stars increase in the 
region of SSC A supporting the hypothesis of an expanding shell. The stellar kinematics closely 
follow the motion of the gas. Analysis of high resolution H I data clearly reveals the presence of 
an H I cloud that appears to be impacting the eastern edge of NGC 1569. Also, an ultra-dense 
H I cloud can be seen extending to the west of the impacting H I cloud. This dense cloud is likely 
the remains of a dense H I bridge that extended through what is now the central starburst area. 
The impacting H I cloud was the catalyst for the starburst, thus turning the dense gas into stars 
over a short timescale, ~ 1 Gyr. We performed a careful study of the spectral energy distribution 
using infrared, optical, and ultraviolet photometry producing a state-of-the-art mass model for 
the stellar disk. This mass modeling shows that stars dominate the gravitational potential in the 
inner 1 kpc. The dynamical mass of NGC 1569, derived from V max , shows that the disk may be 
dark matter deficient in the inner region, although, when compared to the expected virial mass 
determined from halo abundance matching techniques, the dark matter profile seems to agree 
with the observed mass profile at a radius of 2.2 kpc. 

Subject headings: galaxies: individual (NGC 1569) — galaxies: dwarf galaxies — galaxies: starburst — 
galaxies: kinematics — galaxies: spectroscopy 
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1. Introduction 

Dwarf galaxies are the most common type of 
galaxy and are proposed as the building blocks of 
larger galaxies in the early universe. Yet, even the 
basic structure of Magellanic-type dwarf irregular 
(dim) galaxies remains controversial. Studies of 
the distributions of projected minor-to-major axis 
ratios b/a have been interpreted to mean that dim 
galaxies are thick disks with intrinsic ratios (b/a)o 



1.5-3 times that of spirals (|Hodge fc Hitchcock 



1966 



1992 



van den Bergh||l988| |Staveley-Smith et al. 



Yet others have interpreted the distribution 
as evidence that dim galaxies are triaxial ( Binggeli 
fc Popescu||1995 1, only a little less spherical than 
dwarf ellipticals (dE) (Sung et al. 1998 iClearly, 
statistical analyses of axis ratios have not been 
adequate to reveal the true shapes of dim galaxies. 

However, the stellar kinematics of a single 
galaxy, defines its intrinsic shape. The ratio 
^m a x/c z is a measure of how kinematically hot 
a system is and, therefore, reveals its structure. 
Here, Vm ax is the maximum speed of rotation of 
the galaxy and a z is the intrinsic stellar velocity 
dispersion perpendicular to the disk. Elliptical 
and dEs have V m&yi l<J z < 1 (Figure 4-6 in Bin- 



ney & Tremaine||1987 


Pedraz et al.||2002 1 whereas 


spirals have V m3iX /a z > 1, usually 2-5 ( 


Bottema 



We study the stellar and gas kinematics of 
NGC 1569, with the goal of determining the three- 
dimensional shape by measuring V max /(T Z . We also 
aim to understand its recent starburst from an 
analysis of the kinematics. The stellar kinemat- 
ics are derived from long-slit spectroscopy of inte- 
grated light, in the sense that the stars are unre- 
solved, and the high-resolution H I kinematics are 



from The H I Nearby Galaxy Survey (Walter et al. 



2008[ THINGS;). Additionally, we discuss a new 
method of imaging the high-resolution H I data 
cube and apply a double Gaussian decomposition 
of the H I line profiles, a method developed by|Oh| 



et al. (20081, in order to separate ordered motions 
from non-ordered motions. 

NGC 1569 is a nearby post-starburst dim 
galaxy. It has three supermassive star clusters 
(SSCs) that formed at the end of the most recent 
episode of intense star formation that concluded 



approximately 10 million years ago (Hunter et al 
2000 and references therein). This dim galaxy 



has complex gas kinematics. Stil & Israel (19981 
used the Westerbork Synthesis Radio Telescope 
(WSRT) to obtain high-sensitivity H I emission 
maps of NGC 1569 and show the presence of a 



strong non-circular motion (NCM) H I cloud ( Stil 
& Israel 2002 ) and a possible H I companion and 



bridge that are to the east of the disk ( Stil fc Israel 
1998j |Miihle et~aTpM5] ). |Miihle et al.| ( |2005[ ) also 



observed these features and suggest that they are 
being accreted onto NGC 1569. In the center of 
the galaxy, energy from supernova explosions from 
the most massive stars in SSC A has blown a hole 
in the H I gas (Israel & van Driel 19901. There 
is X-ray emission to the north and south of the 
galaxy, perpendicular to the disk, from the outflow 



of hot, ionized gas driven by the starburst (Ott 



et al. 2005a), and strong magnetic fields appear 



to be channeling the gas outflow (Kepley et al. 



2010). There is an Ha chimney extending to the 



north showing the escape route of the gas (Hunter 



et al. 1993 Miihle et al. 20051, which has been 



suggested as evidence for blowout ( |Westmoquette 
etaL][2008l ). 



Figure [T] outlines the major H I features, in- 
cluding the suggested H I companion and bridge 
dStil fc Israel|[2002| , the NCM H I cloud ( |Stil fc 



Israel 2002), an ultra-dense (UD) H I cloud de- 
scribed below, and extended H I emission imaged 
in this work. This figure also identifies the Ha 



chimney and 'western arm' (Hodge 1974), shows 
the location of SSCs A and B, and marks the lo- 
cation, size and position of the optical slit used 
to determine the stellar kinematics along the ma- 
jor axis of NGC 1569 discussed in this work. The 
sense of rotation of the stellar and gas disks are 
marked by the white curved arrows and the near 
and far sides of the disks compared to Earth are 
also labeled. The relative velocities with respect 
to the systemic velocity, V sys , of the H I clouds 
and features are also shown in Figure [T] 

Table [T] gives some global parameters of NGC 
1569. Note that the d istance (3.36 ± 0.20 Mpc; 
Grocholski et al.|2008[ ) places NGC 1569 comfort- 



ably within the IC 342 galaxy group, which has an 
average distance of 3.35 ± 0.09 Mpc ( |Karachent-| 



sev 



2005). Having a group association raises the 
possibility of NGC 1569 interacting or merging 
with other galaxies in the IC 342 galaxy group, 
potentially explaining its unusual characteristics. 
NGC 1569 is, in fact, a complex system with 
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complicated kinematics, and we will present ev- 
idence here and in a forthcoming paper (Johnson, 
in prep) of an interaction or a past merger that 
has created the dim galaxy we observe today. 

This work is part of a large collaborative effort 
called Local Irregulars That Trace Luminosity Ex- 
tremes, The H I Nearby Galaxy Survey (LITTLE 



THINGS; Hunter & Elmegreen|2012 1. One of the 



goals of the LITTLE THINGS survey is to un- 
derstand how dim galaxies form stars. LITTLE 
THINGS is a multi-wavelength study of 37 dim 
and 4 Blue Compact Dwarf (BCD) galaxies. LIT- 
TLE THINGS was granted close to 376 hours in 
2007 and 2008 on the NRACQ Very Large Array 
(VLA) in the B-, C-, and D-array configurations 
in order to obtain deep H I maps for 21 galaxies. 
The other 20 galaxies are from the VLA archives. 
In addition to the H I data, LITTLE THINGS has 
extensive Ha and optical photometry for each dim 
galaxy from previous studies ( |Hunter fc Elmegreenj 



2004 2006) and has also collected GALEX ultra- 



violet images for most and near- and mid-infrared 
data (ground based JHK and Spitzer) for a sub- 
sample. 

Our paper is structured as follows: Section 2 
describes the observations and data reductions for 
the optical spectroscopy and the H I data used 
in this study. Section 3 discusses our data anal- 
ysis; Section 4 outlines our observational results; 
Section 5 describes our mass modeling; Section 6 
is our discussion; Section 7 is a summary of our 
results. 

2. Data 

2.1. Stellar Optical Spectroscopy 

2.1.1. Observations 

We observed NGC 1569 with the Kitt Peak 
National Observatory (KPNCQ Mayall 4-meter 
telescope with the Echelle spectrograph for three 
nights in February 2008 and five nights in Novem- 
ber 2008. The sky conditions for the February 



1 National Radio Astronomy Observatory (NRAO) is a fa- 
cility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 

2 The Kitt Peak National Observatory is operated by the 
National Optical Astronomy Observatory (NOAO), which 
is operated by the Association of Universities for Research 
in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 



2008 observing run were moderate with only a 
few cirrus clouds interfering sporadically through- 
out the run. The moon was in a waning cres- 
cent phase and was separated by more than 180° 
from NGC 1569. The weather conditions for the 
November 2008 observing run were superb with 
nearly photometric conditions throughout the five 
nights, however, the moon passed through third 
quarter during the second night of our run but 
was separated from NGC 1569 by more than 60°. 

Table [2] lists the observation log. In order to 
place the slit, we identified the center of the galaxy 
and position angle (PA) of the major axis from 
iso-intensity contours on our V-band image. We 
observed at four PAs, all centered on the morpho- 
logical center of the galaxy: the major axis (PA = 
121°), the minor axis (PA = 31°) and ±45° from 
the major axis. Figure |2| shows the ffS^V-band, 
O III A5007A, H/3 and Ha filter composite image 
of NGC 1569 taken with the Advanced Camera for 
Surveys and the Wide Field Planetary Camera 2 
instruments (bottom panel) with the slit positions 
superimposed. 

We converted the Echelle into a two-dimensional, 
long-slit spectrograph by replacing the cross- 
disperser with a mirror flat. This permitted the 
isolation of a single order and produced a 3' 
long slit. In order to maximize throughput, while 
maintaining high spectral resolution, we opened 
the slit to a width of 2'.'5. For the stellar kinematic 
data, we targeted the prominent Mg lb stellar ab- 
sorption features at wavelengths 5183.6 A, 5172.7 
A, and 5167.3 A, by isolating order 10 for the 
316-63° Echelle using the post-slit filter KP 1433. 
This filter has a central wavelength of 5204 A and 
a full-width half-maximum (FWHM) of 276 A. 
The free spectral range of the Echelle was 517 A 
centered at 5170 A, but, for the Tektronics 2048 
x 2048 CCD that we used, only 286 A fit across 
the chip. We used the "fast UV" camera along 
with pre-slit filter GG 495. The prominent Mg 
lb stellar absorption features were easily detected. 
There are no telluric features in the spectral region 
we observed. 

For our spectral range, this setup provided a 



3 Hubble Space Telescope (HST), is operated by NASA and 
ESA at the Space Telescope Science Institute (STScI). 
STScI is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 
5-26555. 
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FWHM = 0.38 A as measured from the central 
lines of Th-Ar comparison lamp exposures, giving 
a velocity resolution of ^22 km s _1 . This cor- 
responds to a velocity dispersion, er, of ^9 km 
s _1 . We achieved an angular pixel scale of about 
l'/2 pix -1 , after binning by 2 pixels in the cross- 
dispersion direction, and a final spectral pixel scale 
of 0.14 Apix- 1 . 

For each observing run, we moved a single ra- 
dial velocity standard star in 30" steps along 
the slit and used these data to trace the cross- 
dispersion curvature of the slit on the CCD. We 
averaged 10 dome flats from each night to de- 
termine pixel-to-pixel variations in the detector. 
Also, we took twilight sky exposures that were 
used to correct for illumination variations. The 
rotator on the Echelle was set to the desired PA 
before observations began, and we generally used 
one PA during a night. 

Each night we observed a series of radial veloc- 
ity standard stars. These Milky Way stars were 
selected to have a range of spectral types (F, G, 
and K) and luminosity classes (III - V) to pro- 
vide a varied spectral stellar library for the cross- 
correlation analysis. Table [3] lists the standard 
stars used in the cross-correlation. We observed 
some of the same stars from night to night and 
run to run to minimize and identify systematic 
offsets. There were no trends identified over our 
narrow wavelength region between the radial ve- 
locity standard star spectral types or classes. 

For NGC 1569, we took 1800 s exposures after 
first verifying the telescope pointing by offsetting 
from a nearby foreground star before each expo- 
sure. To ensure proper sky subtraction after each 
galaxy exposure, we observed a nearby dark region 
of sky for 1800 s. In some cases, the sky exposures 
were scaled by a constant to correct for differing 
sky conditions. Following every star, galaxy, or 
sky exposure, we took Th-Ar comparison lamp ex- 
posures in order to acquire accurate wavelength 
calibration as a function of time and position. For 
each PA, we took five to seven exposures on the 
galaxy, which we co-added for a total integration 
time of 2.5 - 3.5 hours, during a single night. 

2.1.2. Data Reductions 

All data reductions were accomplished us- 
ing the Image Reduction and Analysis Facility 



(IRAF£|. We fit the overscan to all images in 
order to remove the electronic pedestal, and cor- 
rected for pixel-to-pixel variations using dome flats 
taken each night (CCDPROC). We trimmed the 
two-dimensional flats so that the steep drop-off of 
light towards the edges would not complicate the 
normalization function. 

The wavelength solution was modeled by fitting 
the comparison lamp exposures independently for 
each spectrum each night. For a given night, the 
wavelength solutions were stable to within a tenth 
of an angstrom. Because we observed comparison 
lamp exposures before every spectrum, our wave- 
length solutions were robust. Thus, we were able 
to place all images on the same wavelength scale 
for each observing run (longslit, identify, rei- 
DENTIFY, fitcoords). We used the radial veloc- 
ity standard star that was moved spatially along 
the slit to map out the distortion along rows of 
constant position, and then, we used the compar- 
ison lamp exposures to trace the distortion along 
columns of constant wavelength. These fitting 
procedures produced coordinate maps that were 
then applied to every respective object exposure 
(transform). Next, we applied a slit function 
that corrected for the minor illumination drop off 
at the edges of the slit in the cross-dispersion 
direction. This was achieved by combining the 
twilight flats from a single night and fitting a 
two-dimensional function to the illumination pat- 
tern along bins of summed columns in the cross- 
dispersion direction (illum). 

Before combining the spectra, we ensured 
proper galaxy alignment by over-plotting the spec- 
tra and matching the peak flux in individual im- 
ages. If necessary, a pixel shift was applied to 
align peaks. Galaxy images from a single night 
with the same PA were combined into a single 
spectrum with an algorithm to reject cosmic rays 
(imcombine). The same combining process, with- 
out the shifting procedure, was applied to the sky 
images. Finally, we subtracted the co-added sky 
image from the co-added galaxy image while tak- 
ing care to scale the combined sky image, when 
needed, to correct for differing exposure times 
and sky conditions. Figure [3] shows the co-added 
galaxy spectrum of the major axis of NGC 1569 



4 IRAF is distributed by the National Optical Astronomy 
Observatory. 
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before and after sky subtraction. 

2.2. Ionized Gas Optical Spectroscopy 

2.2.1. Observations and Data Reductions 

A different post-slit filter (KP f547) was used 
for the ionized gas, but otherwise, the telescope 
setup was unchanged. The KP 1547 filter has a 
central wavelength of about 5015 A and a FWHM 
of 27 A, and was used for targeting the O III A5007 
A emission feature. We observed both the ma- 
jor and minor axes. We took dome flats, twilight 
flats, and stepped one radial velocity standard star 
along the slit spatially to correct for the slit distor- 
tion and to find the center of the slit in the same 
fashion as for the other filter, described previously. 

Additionally, we took Th-Ar comparison lamp 
exposures after all object images for accurate 
wavelength calibration. We observed three galaxy 
and three sky positions for the major axis PA, and 
three galaxy positions for the minor axis. How- 
ever, due to threatening clouds during the course 
of the night, we were only able to observe one 
sky image for the minor axis. The sky image was 
scaled to compensate for the large difference in 
exposure time. Despite the difficulties, the bright- 
ness of the O III emission line made it possible to 
obtain high signal-to-noise emission line spectra. 
We followed the same data reduction procedure as 
for the stellar spectral data. Figure [4] shows the 
fully reduced, two-dimensional, co-added, O III 
A5007 A emission spectrum for the major and 
minor axes of NGC 1569. 

2.3. H I Data 

The LITTLE THINGS project retrieved the 
calibrated H I line Mi>-data from THINGS but 
improved upon the quality of the THINGS data 
products by applying a multi-scale cleaning algo- 
rithm implemented in the AIPS task IMAGR (see, 
e.g.,|Cornwell|2008||Rlch et al.|2008|[Greisen et al.| 
2009[ ). Instead of approximating all structures as 
a sum of point sources, this algorithm iteratively 
fits Gaussians to four angular scales of 135" , 45" , 
15" , and 0" , respectively. The multi-scale clean- 
ing algorithm produces a single data cube in which 
all observed angular frequencies are properly rep- 
resented and with near uniform (Gaussian) noise 
statistics, which is a great advantage over stan- 
dard cleaning methods. A beam size of 7'/71 x 



7" 04 (with a position angle of 16? 5) and a spectral 
resolution of 2.6 km s _1 were achieved. 

In order to separate noise from real galaxy 
emission, we used the Natural- weighted cube and 
smoothed it to a 25" x 25" (convl) beam. We 
blanked this smoothed cube twice, first for all pix- 
els below 2.5 a of the rms of the smoothed cube 
(blank), and then again by hand for any addi- 
tional emission that was clearly noise. This final, 
blanked cube became the master for the cubes 
with smaller beam sizes and was used as a con- 
ditional transfer. For a comprehensive review of 
our mapping procedure, please see Hunter et al. 
(submitted) . 

The half-power beam width (HPB W) of the pri- 
mary beam is 32' , but the largest angular scale 
of emission that is expected to be recovered by 
the interferometer is ^15' , which is determined 
by the D-array. Any smooth structure larger than 
^15' will be invisible to the interferometer. This 
was useful in removing foreground Milky Way HI 
emission, which has velocities that partly overlap 
with some of the velocities encountered in NGC 
1569. In essence, the interferometer is blind for 
emission from the Milky Way, and will only regis- 
ter its small-scale structure, which is still consid- 
erably more extended than that of NGC 1569. The 
multi-scale cleaning algorithm was able to model 
this Milky Way emission sufficiently well so emis- 
sion from NGC 1569 could be isolated from the 
rest using the blanking method described above. 
A more in depth discussion of Milky Way confu- 
sion with NGC 1569 is discussed in our upcoming 
paper describing the map made with the Green 
Bank Telescope (GBT) (Johnson, in preparation, 
hereafter Paper II). 

Figure [5] shows the LITTLE THINGS inte- 
grated flux map. This image shows tenuous emis- 
sion to the south and northeast of the galaxy that 
are imaged for the first time. Although faint, 
this low-level, extended emission was confirmed 
in follow up observations with the Robert C. Byrd 
Green Bank Telescope (GBT) and will be dis- 
cussed in Paper II. 
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3. Analysis of Data 

3.1. Stellar Optical Spectroscopy 

3.1.1. Extracting One- dimensional Spectra 

We extracted one-dimensional stellar spectra 
from the two-dimensional spectra of NGC 1569 
(scopy). It was necessary to sum over a series of 
consecutive rows in the cross-dispersion direction 
during spectral extraction in order to increase the 
signal-to-noise (S/N). In the portion of the spec- 
trum that contains the brightest light, from radii 
0" to -25" (West of center) measured from the 
center of the galaxy, which is also where the SSCs 
reside (shown in Figure [6| , we summed over five 
rows, 6'.'1 of the major axis. For the regions of the 
galaxy that have a lower light concentration, from 
radii 17" to 0" (East of center), we summed over 
ten rows of the major axis, 12'.'2. To avoid diluting 
the true stellar velocities and velocity dispersions, 
we never summed over more than 12'.'2, which, at 



a distance of 3.36 Mpc (Grocholski et al. 20081, 
corresponds to 195 pc. 

Figure [7] shows three, one-dimensional galaxy 
spectra and one stellar template spectrum for com- 
parison. The top and third panels were the result 
of summing over 12" 2 in the cross-dispersion di- 
rection along the major axis of NGC 1569 and the 
second panel summed over 6'/l. These three spec- 
tra show the range of S/N for our galaxy spectra 
and all of our spectra used in the CCM had a S /N 
above ~ 5. The bottom panel shows HD 132737, 
a K0 III template star for comparison. All of the 
standard stars had S/N above ~ 50, usually near 
90, although this was difficult to measure due to 
the lack of line-free continuum (see bottom panel 
of Figure ^ . 

We extracted one-dimensional spectra of the ra- 
dial velocity standard stars using APALL (see, e.g., 
bottom panel of Figure [7]) . We identified the cen- 
ter of the peak emission in each star and fit the 
background using the spatial regions of the slit 
that contain no stellar data. In some cases, we 
used the twilight sky exposures to extract a solar 
spectrum, which was used in the cross-correlation 
method (CCM) as a G2 V radial velocity standard 
star. Table [2] identifies the nights that utilized 
these solar spectra by the +1 added to the number 
of standard star spectra in column 9. Once one- 
dimensional spectra were extracted, we added into 



the header of each star the Vheiio radial velocities 
from the United States Nautical Almanac (2008). 
These spectra of radial velocity standard stars are 
the templates that were used in the CCM. We fit 
a cubic spline function to the continuum in the 
spectra of the galaxy and the template stars (con- 
tinuum). The slit width of 2'.' 5 may have caused 
small telescope tracking errors for some of the tem- 
plate stars because the seeing on most nights was 
around 1"2. However, the exposure times for the 
standard stars were short, less than a minute for 
a large portion (see Table |3j, so tracking errors do 
not a play a large role in the total error budget. 

3.1.2. Cross-Correlation Method 

We cross-correlated the template stars with 
each other to test the CCM and to check the ob- 
served stars' velocities, which was within a few 
percent of the stars' known velocities in most cases 
(fxcor). The one-dimensional spectra from the 
galaxy were cross-correlated with the stellar tem- 
plates observed during the same night. We only 
performed cross-correlations of the stars with stars 
and stars with galaxy spectra that were obtained 
on the same night in order to minimize the uncer- 
tainties in the observational setup, such as varying 
position angles. The heliocentric radial velocities 
and velocity dispersions of the stars in NGC 1569 
were determined, for each spectrum, by fitting a 
single Gaussian curve to the peak of the CCM. 
We did not explore the possibility of asymmetric 
Gaussian profiles as the quality of our data were 
not good enough for accurate analysis. We made 
the assumption that the stars in NGC 1569 are 
virialized and the Gaussian profiles are symmet- 
ric. We excluded regions around emission line fea- 
tures in order to avoid internal galaxy emission or 
residuals due to atmospheric emission lines that 
were not completely removed during the sky sub- 
traction step (see the top image of Figure [3] which 
shows residual sky emission lines around 5196.8 A 
and 5199.1 A, and the top panel of Figure [7]). 

We cross-correlated each template star with 
each galaxy position and averaged the Vheiio val- 
ues by weighting them by the standard deviations 
of the different template spectra, thus determin- 
ing the Vheiio values at each galaxy radius. By 
weighting the Vheiio values in this manner, an ill- 
fitting template would not contribute as much to 
the resulting Vheiio of the galaxy spectra MS (\ good- 
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fitting template. Also, we observed no trend be- 
tween stellar spectral types or classes when deter- 
mining a "good-fitting" versus an "ill-fitting" tem- 
plate in the CCM; S/N dominated the quality of 
the fit. 

To determine the intrinsic velocity dispersions 
of the stars, <7 bs, in NGC 1569 as a function of ra- 
dius, we employ the definition of a Gaussian and 
correct for instrumental broadening using the fol- 
lowing equation: 

v/FWHM obs 2 - FWHM instr 2 
CTobs - 2^5 (1) 

We determine the inherent, average instrumental 
full-width half-maximum, FWHMi nstr = 35 km 
s _1 , by cross-correlating the template stars from a 
single night with each other. The intrinsic broad- 
ening of a single star should be unresolved, so the 
FWHM instr of the cross-correlation of the tem- 
plate stars with each other represents the instru- 
mental signature. Although the spectral types and 
classes of the radial velocity standard stars span 
a range, our usable wavelength region is only ~ 
150 A, thus cross-correlations of the stars between 
differing spectral types did not create variations 
in the resulting FWHM instr or FWHM obs above 
their respective standard deviations. 

3.2. Ionized Gas Optical Spectroscopy 

3.2.1. Extracting One-Dimensional Spectra 

Due to the starburst nature of NGC 1569, many 
star forming regions lie throughout the observable 
disk and create bright O III features seen along the 
cross-dispersion direction of these spectra. Figure 
[4] shows the two-dimensional spectra of the O III 
A5007 A emission feature observed in NGC 1569 
along the major axis, top panel, and minor axis, 
bottom panel. Because the emission features are 
so bright, we are able to extract a single spectrum 
from each cross-dispersion row, giving an angu- 
lar resolution of 1'.'09 per pixel (17.7 pc). Also, 
we achieve a velocity resolution of ^27 km s _1 as 
measured from a comparison lamp emission line 
at A = 5014 A with a FWHM of 0.46 A. This cor- 
responds to a velocity dispersion, a, of ~12 km 



3.2.2. Multiple Gaussian Decomposition 

Figure [8] shows an example of a single extracted 
spectrum of the O III A5007 A emission feature 
taken from the major axis. It is evident that there 
are three distinct, bright peaks over a total width 
of 4.1 A, or, ~176 km s _1 . We disentangled these 
blended features with multiple Gaussian compo- 
nents (splot). The three velocities obtained for 
the brightest three peaks in Figure [8] are -133.9, 
-84.2, and -31.5 km s _1 spanning a range of over 
100 km s _1 at a distance of -13'.'2 from the center of 
the galaxy along the major axis. This fitting pro- 
cedure was repeated for each spectrum extracted 
along the cross-dispersion axis. 

3.3. Extracting the H I Bulk Velocity 
Field 

3.3.1. Tilted Ring Models 

We use the newly imaged LITTLE THINGS H I 
data to determine the maximum rotation, y max , of 
the H I gas for NGC 1569. One common method of 
constructing rotation curves of disk galaxies from 
high-resolution H I data is the tilted ring model 
(Begeman 1989). The first step in applying this 
model is to create a two-dimensional velocity field 
out of a three-dimensional H I data cube. Then, 
from the extracted velocity map, a series of ring 
parameters are determined by independently fit- 
ting a series of concentric, circular rings in the 
plane of the disk. The ring parameters are: in- 
clination of the galaxy disk (z); PA of the ma- 
jor axis (9) measured east from north; kinematic 
center (X pos , Y pos ); systemic velocity (V sys ); ex- 
pansion velocity in the radial direction of the disk 
(K.xp); and, finally, the rotational velocity (V ro t)- 
The ring size and spacing between rings must be 
fixed and are usually determined using the beam 
size as the smallest ring radius and ring width. In 
the tilted ring model, these rotational parameters 
are determined by assuming an axisymmctric thin 
disk with closed circular orbits and are fit inde- 
pendently to each of the pre-determined rings in 
succession. 

To determine the rotation curve of NGC 1569 
we use the GlPS^jtask rotcur, which applies 
the tilted ring model. We input the bulk velocity 



5 The Groningen Image Processing SYstem (GIPSY) has 
been developed by the Kapteyn Astronomical Institute. 
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field, described in the following section, along with 
the rotation parameters derived after eight itera- 
tions, which is how many iterations were required 
for the model to converge. We used rings that were 
8" wide (roughly one beam width) and separated 
by 8" so that there was no overlap between rings. 
Our final ring radius was 200" because beyond 
this radius the signal-to-noise dropped drastically. 
The final parameters determined from the tilted 
ring model are as follows: the kinematic center, 
(A pos = 04:30:46.125, y pos = +64:51:10.25), which 
is marked by the blue star on Figure [2] and differs 
by 58" from the morphological center used for the 
stellar kinematics; i = 69 ± 7; 8 = 122°, which is 
only one degree different from the PA of the major 
axis used for the stars; V sys = -85 km s _1 , which 
agrees within the error of the systemic velocity de- 
rived for the stars; V CKp = km s" 1 . Table [T] lists 
these final parameter values determined from the 
tilted ring model. All velocities for the H I data 
are heliocentric velocities, the same velocity refer- 
ence frame as the stellar data. 

3.3.2. The Bulk Velocity Field 

NGC 1569 shows a dominance of non-circular 
and random motions as seen in the intensity- 
weighted mean velocity field in Figure [9] The blue 
features across NGC 1569 show H I gas with mo- 
tions that are non-circular and do not follow the 
bulk motion of the galaxy disk. A 'typical' rotat- 
ing disk generally displays a grading scale of isove- 
locity contours with blue-shifted (relative to the 
systemic motion) H I on one side and red-shifted 
H I on the opposite side. The systemic velocity of 
the system separates the two halves. This is not 
what is observed in Figure [9] Nevertheless, we 



were able to apply the method developed by Oh 



et al. (2008) to extract bulk rotation from chaotic 



systems like NGC 1569. 

An iterative double Gaussian deconvolution 
procedure was used to separate velocities of H I 
clouds superimposed upon the bulk motion at a 
given position in the H I data cube. The deconvo- 
lution program uses only the H I intensities that 
are above 2 a in the bulk velocity field. This 
limits the extracted velocity fields to the inner, 
concentrated H I disk. The red outline in Fig- 
ure [9] roughly identifies the region over which we 
extract the bulk velocity field. At each position 
within this region, if there are two distinct Gaus- 



sian peaks detected above 2 a, then the velocities 
of the central wavelengths of each respective peak 
are compared to a model velocity field. The veloc- 
ity closest to the model velocity field is assigned to 
a bulk motion map while the outlying velocity is 
assigned to one of two non-circular motion maps. 
The strong non-circular motion map contains all 
velocities that are derived from Gaussian peaks 
that have a higher intensity than the bulk motion 
velocity, whereas the reverse is true for the weak 
non-circular motion map. 



Figure 10 shows (a) the final bulk velocity map, 
(b) the strong non-circular velocity map, (c) the 
weak non-circular velocity map, and (d) the strong 
non-circular motions superimposed on the weak 
non-circular motions map. The black oval in each 



panel of Figure 10 identifies an H I cloud that 
has systemic motions significantly deviating from 
the bulk motion of the galaxy. In fact , this cloud 
is virtually absent in the bulk velocity field (see 
panel (a) in Figure |j"o|. |Stil fc Israel] ( |2002| ) first 



detected this H I cloud and describe it as "high 
velocity cloud 3 (HV3)". This cloud is discussed in 
more detail in Section [5] In Section [4. 3| we derive 
the H I rotation curve using the two-dimensional 
bulk velocity field and the tilted ring model. 

4. Results 

We cannot assume that the stellar and H I disks 
have the same geometry because there are exam- 
ples of dim systems in which the stars and gas 
have different kinematics. In the case of WLM, 



Leaman et al. ( 2009 1 compare the kinematics of 



the stars and H I and discover that the stars and 
H I are kinematically decoupled and do not pos- 
sess the same geometry. They conclude that the 
stars lie in a thicker disk than the H I. IHunter et all 



(|2002j) determined that the stars in NGC 4449 are 
in a face-on disk, whereas the H I disk is inclined 
and that it exhibits rotation as well as a counter- 
rotation component in the center. Therefore, we 
are careful here to examine the stellar kinematics 
and gas kinematics independently. We begin by 
investigating the stars alone. 

4.1. Stellar Optical Spectroscopy 

Figure 



11 



shows position- velocity (P-V) graphs, 
where the positive positions are in the eastern di- 
rection. The error bar in the radial direction rep- 
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resents the angular area of the slit over which the 
velocities were determined. We used the weighted 
average of the minor axis velocity points to deter- 
mine a systemic velocity, V sys , of -82 ± 7 km s _1 , 
which is denoted by the horizontal solid line in 
Figure [TT] This va l ue is consistent with H I stud- 
ies such as Reakes ( 1980 I, V S y S = -77 ± 1 km s _1 
and de Vaucouleurs et al. (1991) 
km s 



are well matched (see Section 4.2 1, from here on 



V s 



sys 



-89 ± 5 



But, it is inconsistent with the value ob- 
tained by |Schneider et al.| ( fl992] ), V sys = -104 ± 4 
km s _1 , who state that the H I flux used to derive 
the velocity is confused with foreground Galactic 
emission. 

We used the four PAs to constrain the inclina- 
tion and line of nodes. The de-projection was de- 
rived by the following formula obtained from pri- 
vate communications with Vera Rubin: 



V S y S — 



V(R) sin(i) cos(t? - 9) 
sec 2 (i) — tan 2 (i) cos 2 (ry — 9)] 1 / 2 



(2) 

Here, i is the inclination of the disk, r\ is the PA 
of the observed point in the plane of the sky, 9 is 
measured north through east to the receding ma- 
jor axis, which is the PA of the major axis (line of 
nodes) , V^bs is the observed line-of-sight radial ve- 
locity, Vgys is the systemic velocity, and V(R) is the 
rotational velocity. Although the velocities of the 
stars are quite irregular, we were still able to con- 
strain the inclination of the disk to i ~ 60° ± 10° 
and the major axis PA « 120° ± 8° using a least 
squares fit. Both of these estimates agree well with 
the H I results (see Table [lj. Furthermore, Fig- 
ure [12] shows P-V diagrams of the stars plotted 
with H I intensity contours. These plots demon- 
strate the kinematic similarity between the stars 
and H I, as the stars follow the H I in all four PAs. 
The stellar velocity dispersions are shown in Fig- 
ure [13] for each of the four PAs. The average of 
the velocity dispersion, < er b s >, weighted by the 
uncertainty of each point, for all four PAs is 21 ± 
4 km s _1 . 

In order to disentangle <tr, <t z , and from 
Cobsj we use the following relation from |Gerssen| 
efaT|p97|: 



o'obs = [°r sin2 ( ? 7-^)+^ cos 2 (77-6*)]sin 2 i+cr z cos 2 i 

(3) 

Here, 77 is the intrinsic angle within the plane of 
the galaxy, 9 is the receding major axis PA and i 
is the inclination of disk. Since the H I and stars 



we use the better determined inclination of i = 
69° ± 7° and major axis PA of 122? 5 from the H I 
kinematics. Since we have only a few points along 
each intermediate PA, i.e., ± 45° from the major 
axis, we focus on the major axis, where 77 = 9 and 
thus, or does not contribute to cr b s . Following 
Swaters (1999), we make the assumption that the 
stellar velocity dispersion in the stars follows the 
Milky Way properties observed for stars near the 
Sun. Therefore, <7 b s = cr z for the major axis. 

Swaters (1999) found a relation between inte- 
grated Mb of a galaxy and the central vertical 
velocity dispersion: 



log a z 



-0.15A/ B - 1-27 



(4) 



This relation predicts a a z of 21 km s" 1 for NGC 
1569, which is equal to our observed value of 21 ± 
4 km s -1 . NGC 1569 agrees well with the overall 
trend observed in spiral galaxies, i.e., lower lumi- 
nosity systems have lower velocity dispersions. 

4.2. Ionized Gas Spectroscopy 

Figure [14] shows the decomposed line-of-sight 
radial velocities of the ionized gas for the major 
axis and has the stellar and H I velocities over 
plotted for comparison. The H I velocities in this 
figure come from a P-V slice through the major 
axis of the H I data cube with a sampling width of 
8" , equal to the FWHM of the H I beam. There 
was no deconvolution done to produce these H I 
velocities; they are simply the observed velocities 
along the major axis position angle. The velocities 
of the O III ionized gas match the velocities of Ha 
along the major axis observed by |Tomita et al.| 
(1994 see their Figure 3). 

In Figure \14\ we observe two forms of ionized 
gas: high velocity gas and low velocity gas. The 
low velocity ionized gas has velocities that match 
the H I and stellar velocities and is likely asso- 
ciated with HII regions, or the remnants thereof, 
leftover from recent star formation. The high ve- 
locity ionized gas at a radius of -20" has the signa- 
ture of an expanding shell as identified by the solid 
red outline in Figure [14] This expanding shell has 
an expansion velocity of ^70 km s _1 . The stellar 
velocities over plotted in Figure[l4]perfectly match 
the velocities of the low velocity gas (~ +10 km 
s _1 ) at the center of the expanding gas shell (R ~ 
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—20" ). These low velocities of the ionized gas also 
coincide with the bulk motion of the H I gas. Cu- 
riously, this agreement extends to the peculiar dip 
to lower velocities at the western edge of the red 
oval outline. This agreement in velocity (stars and 
O III) shown by the black arrow in Figure [l4j are 
presumably from HII regions connected with par- 
ticularly bright stars. Within this expanding shell, 
there is only one H I velocity point. The absence 
of H I in this region is due to an H I hole blown 
out by SSC A QIsrael fc van Driel|p90| . SSCs A 
and B reside at radii -22" and -12" , respectively. 
Due to the lack of H I and the matching stellar 
and low velocity ionized gas in the expanding shell 
region, even including peculiar motions, we con- 
clude that our stellar spectra are dominated by a 
young, luminous stellar population that are ioniz- 
ing nearby gas left over from the star formation 
process. At the same time these young, luminous 
stars, mostly from SSC A, created the expanding 
shell from winds and supernova explosions during 
the most recent star formation episode. This ex- 
panding shell was also observed and discussed by 
Westmoquette et al. ( 2007a|b l, who also observed 
both the low and high velocity ionized gas. 

The high velocity ionized gas seen everywhere 
else is likely the result from outflow from the recent 
starburst activity (Westmoquette et al. 2007a|b I. 
We also see evidence of galactic wind at radii 
greater than +0" . Because of the exponential 
drop in stellar light, the integrated emission from 
the stars was too faint to extract stellar kinematic 
information in our optical spectra beyond a radius 
of ±30" . Thus, no stellar velocity information is 
known in the region of the high velocity ionized 
gas at radii greater than +30" . 

The low velocity ionized gas centered at R = 
+15" in Figure 14 also shows matching stellar 



and gas velocities. This region, however, has more 
H I than the expanding shell centered at —20" . 
It is not as obvious here as for the expanding shell 
centered at R = —20" , but it is possible that the 
stellar spectra around R = +15" are also domi- 
nated by young, luminous stars because they, too, 
follow the velocities of the low velocity ionized gas 
and H I. 



Studies of HII regions of NGC 1569 ( 


Greve et al. 


1996 


Hunter et al.||1993| Waller 1991; Gotz 1990 


Young et al.||1984| show that there is 


star forma- 



intense as the western side of the galaxy where the 
SSCs reside. Again, we note that the velocities of 
the ionized HII regions (low velocity gas) and H I 
agree well in this area. 

4.3. H I rotation curve 

Figure [15] shows the observed rotation curve. 
We find a maximum rotational velocity, Vf t, of 
33 ± 10 km s" 1 , which was calculated from the 
average of the last nine data points in Figure 15 
The uncertainty is dominated by the rapid drop 
in signal-to-noise in the bulk velocity field beyond 
a radius of 2.2 kpc. Within the errors, Figure 
[15] shows a rising, solid body-type rotation curve, 
which at best reaches the usual flat part at radii 
> 2.2 kpc. 

Figure[T6]shows a velocity dispersion map of the 
bulk (main) velocity component only. The average 
H I velocity dispersion is high, < ohi > = 10 
km s j compared to the observed rotation speed, 
V TO t = 33+10 km s _1 . Because of the high H I 
velocity dispersions in NGC 1569, we correct the 
observed rotation velocities for asymmetric drift in 



the same manner as Oh et al. (20111, and derive 



tion in the region at R = +50" , but it is not as 



the curve shown by the filled circles in Figure [17] 
The maximum rotation velocity determined from 
H I, increased from the observed V TO t = 33 km s _1 
t° Vmax = 50 km s _1 after the asymmetric drift 
correction was applied. 

5. Discussion 

5.1. The Inner Kinematic Structure of 
NGC 1569 

The upper left panel in Figure [12] displays the 
observed radial velocity of gas (contours) and in- 
tegrated stars (points) function of position 
along the slit that passes through the major axis. 
This P-V diagram shows "forbidden" H I motion in 
the upper right quadrant between -40 and -70 km 
s _1 at around a distance of -0'6 from the galaxy 
center, as identified by the black oval. We define 
"forbidden" H I as motion that does not follow the 
ordered rotation. Along the major axis in Figure 
[T2] we show a black curve that traces the aver- 
age velocities of the H I contours from 0' to 2' , 
which is the region of the diagram that follows the 
ordered motion of the H I gas. This black curve 
is then reflected through the origin to show what 
the expected ordered motion should be from -2' to 
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0' . At some positions between distances of 0' and 
-0^5 from the galaxy center, the stellar velocities 
follow the velocity trend of the forbidden motion, 
whereas at other positions around -0.'4, the stellar 
velocities follow the ordered rotation of the H I 
disk. 

Could the gas and stars be showing evidence 
of counter-rotation? If this were the case, then 
we would expect H I contours in the lower left 
quadrant of all the panels of Figure [12] which we 
do not observe. Therefore, we identify the for- 
bidden H I motion as the northeastern tip of the 
NCM cloud identified in Figure [T] and by the black 
oval in Figure [TO] The major axis, PA = 122° 
(shown in Figure 111 and by the black line through- 
out Figure 10), slices through the tip of the NCM 
cloud. The velocities of this cloud in the strong 



ure 13 the velocity dispersion is larger by ~10 km 



non-circular motion map, panel (b) of Figure 10 



match the velocities of the contours of the forbid- 
den H I in Figure [12} The NCM cloud has similar 
velocity structure and position as HV3 in Stil & 
Israel (20021. We conclude that the stars and gas 



are not counter-rotating but, there is evidence of 
a kinematic disturbance. 

When the stellar velocities are laid on top of 
an optical image, we can physically see the areas 
in the galaxy that correspond to observed stellar 
kinematic features. Figure [6]shows the HST image 
in greyscale. Red and blue ovals are superimposed 
on the picture, representing the positions of the ex- 
tracted spectra along the major axis. The length 
of each oval represents the angular area of the slit 
that we summed, while the width of the ovals rep- 
resents the 2'.'5 slit width. Red and blue ovals 
represent red- and blue-shifted Vhdio radial veloc- 
ities, respectively, relative to the systemic velocity 
of NGC 1569. Each of the ovals are representative 
of independent data points from spectra taken on 
different nights. Therefore, there are some red- 
and blue-shifted velocities shown at the same po- 
sitions. The green 'X' marks the center of the slit 
and the center of the galaxy. 

In the region west of SSC A, the overall rota- 
tion of the stars is actually red-shifted, which is 
seen from the dominance of red ovals in Figure 
[6] The locations of the blue ovals that reside to 
the west of SSC A correspond to a distance of - 
27" from the galaxy center. If this location is 
compared with the velocity dispersion plot of the 
stars for the major axis PA (top left panel) in Fig- 



s from the average stellar velocity dispersion. 
The increase in stellar velocity dispersion in this 
region is consistent with secondary star formation 
that formed out of the surrounding material that 
was shocked by the expanding shell discussed in 
Section |4.2| Therefore, these peculiar kinematics 
observed in the stars are likely from light com- 
posed mostly of young, luminous stars in NGC 
1569. 

In the optical HST image of NGC 1569, shown 
in the right image of Figure [I] we can observe Ha 
streamers, which are material likely trapped in the 
outflow, extending to the northeast and southwest 
from the SSC region, as indicated by the arrows. 
Also, there is a large Ha hole where the SSCs re- 
side, also identified in Figure [l] We estimate the 
projected surface area of this hole to be about 0.2 
kpc 2 at the distance of NGC 1569. Figure 12 also 
shows a hole in the H I in all four PAs. It is cen- 
tered around 0'A in the major axis PA and around 



Ofl in the remaining three PAs in Figure 12 (see 
Israel fc van Driel|1990 1 . It is in this hole that we 
observe the disturbed stellar kinematics. These 
features are likely evidence of outflow from me- 
chanical energy input from stellar winds and su- 
pernova remnants from the collection of massive 
stars in the SSCs. Large amounts of X-ray emis- 
sion coming from this region ( Ott et al.|2005a I are 
consistent with this interpretation. 

5.2. Shape of NGC 1569 

We explore the three-dimensional shape of 
NGC 1569 using our kinematic measure, V max /cr z . 
For this relationship, we use the asymmetric drift 
corrected y max from our H I analysis. Because we 
want the maximum velocity of the baryons that 
are responding to the total gravitational potential 
in the galaxy, it is necessary to correct for radial 
pressures in the gas that can artificially mask the 
true maximum rotation speed. We derived an 
asymmetric drift corrected V max of 50 ± 10 km 
s -1 from the H I analysis, discussed in Section 



4.3 and we use cr„ of 21 ± 4 km s from the 



stellar kinematics, discussed in Section 14. 11 We 



follow the same reasoning as Hunter et al. (20051 
for our decision to use <r z in our V/a relation. 
We determine a V max /cr z = 2.4 ± 0.7. This is 
indicative of a thick distribution of young stars 
that dominate our spectra. As a reminder, this 
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result assumes that the stars are in a virailized 
disk and that the Gaussian profiles are symmet- 
ric. Our stellar template spectra span a range of 
classes and types from F-type main sequence to 
K-type supergiants and thus, cover a range of stel- 
lar ages. Our observed galaxy spectra are a mix of 
all the stars in the galaxy but, are dominated by 
the brightest stars, which are the younger stellar 
population. There may be a thin distribution of 
older stars that our spectra simply were not sen- 
sitive to and perhaps, the disk is thicker in the 
region of SSC A than in the center of the galaxy. 
The increase in velocity dispersion and decrease in 
rotation speed around SSC A are consistent with 
velocities of stars that formed from shocked gas 
surrounding the observed expanding shell in the 
region of SSC A, as discussed above. Overall, the 
disk of NGC 1569 appears to be thick compared 
to spiral galaxies. 

5.3. The Outer Morphology of NGC 1569 

Here, we compare the H I morphology and ve- 
locity fields with our deep Ha image of NGC 1569 
and the HST image. The left image of Figure [T] 
shows the contoured integrated H I map. This 
image displays tenuous emission, detected for the 
first time, extending to the south of the galaxy 
disk in a teardrop shape, and to the northeast; 
black ovals outline this emission. Also identified 
in this image are two features to the east labeled 
as 'H I Companion' and 'H I Bridge' identified by 

strife 



Stil & Israel 


(1998) and 


Miihle et al. 


(2005) 



Israel conclude that these features are indicative 
of a low mass H I companion to NGC 1569. Our 
data detect this H I cloud, which can also be easily 
seen in Figures [5] and [9] at around a = 04:32:15, 
5 = +64:49:00. The extended H I emission, cir- 
cle d in Figure |T] and the two features identified 
by Stil & Israel ( 1998 1 as evidence of a possible 



low mass H I companion, may be related to each 
other and show evidence of a recent interaction or 
merger with NGC 1569. We discuss these claims 
further in Paper II where we present results from 
deep GBT H I mapping around NGC 1569 and 
three of its nearest neighbors. 

Figure [18] shows our deep Ha image with the 
contoured integrated H I map overlaid so that we 
can compare the morphologies. What is interest- 
ing here is that there is an UD H I cloud to the 
west that is at the edge of the optical galaxy in the 



same region where we observe the disturbed stellar 
kinematics. This H I cloud is nearly 50 M Q pc~ 2 
more dense than the surrounding medium and if 
there is molecular gas, the column density could 
be much larger. 



In Figure 18 we also observe two smaller, dense 
H I clouds on the opposite side, east of the optical 
galaxy. Possibly, there was one continuous, dense 
H I structure that connected these two dense H I 
clouds with the UD H I cloud. The most concen- 
trated stellar light is in between these dense H I 
clouds where the SSCs reside, which suggests that 
the stars in this region formed out of the dense 
H I ridge that used to connect these clouds. The 
trigger for the starburst could have come from the 
impact of the in-falling NCM H I cloud, shown in 



Figure 10 and identified by the blue oval in Fig- 
ure |18| that compressed the interstellar medium 
of NGC 1569 and caused the dense H I ridge to 
collapse. It is more likely that the NCM cloud 
is in-falling into the galaxy, rather than outflow- 
ing, because there is evidence of shocked gas near 
the impact area of the NCM cloud as seen in the 



continuum maps by |Kepley et al. (2010 see their 



Figure 12). 

We see in Figure [T8| that , in general, the H I 
density is more concentrated to the north of the 
SSCs than it is to the south, as shown by the blue 
arrows and contour levels. There is a slight H I 
depression in the contour levels right above SSC 
A as indicated in Figure [18] possibly created by 
the outflow of the gas. This figure also shows the 
approximate location and projected size of two gi- 



ant molecular clouds (GMCs) identified by Taylor 
et~aT] (fl999"l see Figures 3 and 4). These GMCs 



lie between the UD H I cloud and the SSCs and 
are at the tip of the NCM H I cloud. The region 
between the UD H I cloud and the SSCs that con- 
tains the GMCs is the most probable area for the 
next generation of star formation. McQuinn et al.| 
(2010) suggest that the starburst episode in NGC 



1569 is ongoing and started roughly 500 Myr ago. 



Also in Figure 18 are streamer-like, arced Ha 
features, which trace HII regions, toward the south 
and north. The high density H I in the north 
and the visible Ha features in the south could 
be the results of the projection effects from the 
inclination of the galaxy. The X-ray emission 
is less absorbed to the south, according to |Ott| 
et al. (2005a I, indicating that the north side of 
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the galaxy is pointed toward us, thus closer than 
the south side, assuming that the emission is es- 
caping in the z-direction. Ha arcs and streamers 
extend nearly all the way through the H I disk in 
the southern direction of the disk. Similarly, the 
H I gas appears more dense in the northern half of 
the disk because there is more gas along the line- 
of-sight than there is in the southern half of NGC 
1569. 

Figure [18] also shows a strong Ha chimney to 



the north (Hunter et al. 1993 Miihle et al. 2005 



Westmoquette et al. 2008). The H I depression 



noted on Figure 18 shows the escape route for 



the Ha emission ( |Muhle et al. 2005 Westmoquette 
et al.||2008[). Magnetic field lines found by Kepley 



et al. (20101 may be following the outflow of the 
gas in this region and are channeling cosmic ray 
electrons out of the disk. Although the H I is more 
concentrated to the north, due to the projection 
effect of the inclined disk, the intense stellar winds 
from the star burst are able to penetrate all the way 
through the H I disk in the northern half of the 
galaxy, too. Tenuous Ha emission to the northeast 
of the optical galaxy stretches all the way through 
the H I disk just as it does in the southern part of 
the galaxy. 

Westmoquette et al. (2008) discuss the kine- 



matics of the Ha emission and suggest that 
'blowout,' defined here as gas escaping the gravita- 
tional potential of the galaxy, from stellar feedback 
from the massive stars created in the starburst, is 
not yet occurring due to decreasing ionized gas 
velocities at the edges of the expanding Ha shells 
seen near SSC A. If we assume the total mass 
of NGC 1569 is equal to our dynamical mass of 
1 x 10 9 M (discussed in Section 6.3 1, then the 
escape speed is ~70 km s" 1 at a radius of 2.2 
kpc. According to Figure 14 the ionized gas we 
observe with our O III spectra show high velocity 
gas in excess of 70 km s _1 indicating that perhaps 
some of the gas is in the process of being expelled 
from the galaxy in a blowout 



Ott et al. (2005b I 



also suggest that NGC 1569 is in the process of 
losing about 5 x 10 6 M© of hot gas, corrected for 
our distance. However, Ott et al. argue that in 
order for gas to escape completely from a galaxy, 
it must be moving at supersonic velocities when it 
hits the cool surrounding gas at one scale length. 
We observe tenuous H I stretching nearly 12 kpc 
south of the galaxy that is likely associated with 



tidal filaments (Paper II), so, the hot gas may cool 
and fall back onto the disk of NGC 1569 when it 
hits this cold, neutral gas. Thus, it is not clear if 
this hot gas will actually escape the galaxy. 

Figure [19] shows the H I intensity- weighted ve- 
locity contours, with the extended H I trimmed, 
superimposed over a low contrast Ha image. In 
the region around the Ha image, we observe a but- 
terfly shape in the H I. The western wing of the 
butterfly has velocities that are blue-shifted rela- 
tive to the systemic velocity, as seen by the blue 
and purple contours, whereas the eastern butterfly 
wing has red-shifted velocities, as seen by the or- 
ange and red contours and are identified in Figure 
19 The location of the UD cloud (solid black oval) 
and NCM cloud from Figure [To] (dashed black 
oval) are also indicated. The velocity structure 
of the NCM cloud appears different in Figure [19] 
than in Figure [10] because we are observing the 
intensity- weighted velocity field instead of the de- 
convolved velocity structure. The possible 'H I 
Companion' and 'H I Bridge' discussed by |Stil fc 



Israel ( 1998 1 are seen as purple contours and do 
not follow the rotation of the rest of the galaxy. 
Again, we surmise that these features are most 
likely remnants from a past interaction. 

To the south of the Ha image in Figure [19] 
we see a cavity in the H I. This cavity is likely 
carved out by ultraviolet photons, winds, and ex- 
plosions from the massive stars in the starburst 
X-ray emission fills the cavity ( Martin et al.j 2002 



Ott et al. 2005a I. Furthermore, Figure 19 shows 
that the spectacular Ha tail, known as the 'west- 
ern arm' from Hodge (1974), that extends to the 
southwest of the optical galaxy is, in fact, the lit up 
edge of the UD H I cloud. This western arm was 
previously believed to have been associated with 
the strong NCM cloud created from outflow due 
to the apparent spatial overlap of the two features 



(Waller 


1991 


rael 


( 


2002 1 ar 



However, Stil & Is- 



infall rather than outflow and that these struc- 
tures are connected to the H I bridge. Similarly, 



Miihle et al. ( 2005 1 argue that H I gas is being ac- 



creted onto NGC 1569 and |Ott et al] ( |2005a| | also 
see the Ha tail emitting soft X-rays. Mailing H I 
gas is a likely scenario, especially since we can see 
tenuous H I emission to the south and northeast 
of NGC 1569, which has the same radial veloci- 
ties as the H I bridge. Because the tenuous H I 
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gas to the south is spread over such a large region 
of space, ~92 kpc 2 , and evidence of interaction is 
observed in our GBT data, which we report in an 
upcoming paper, the origin of this tenuous gas, H I 
companion, H I bridge and Ha tail coming from 
outflow from the starburst is unlikely. We con- 
clude that these features are the result of a recent 
interaction, which caused the NCM cloud to infall 
and compress a dense ridge that stretched across 
the entire optical disk, creating a global starburst 
which created the SSCs, and left behind what is 
now the UD H I cloud. 

We see that the velocities of the UD cloud are 
fairly constant, but, according to Figure |16[ the 
velocity dispersions at this location are high com- 
pared to the average. Thus, we observe that the 
UD H I cloud is kinematically warm because of its 
high velocity dispersion. Perhaps this dense H I 
region to the west of the optical disk will be the lo- 
cation of a future star formation episode, provided 
the gas kinematically cools over time. 

6. Mass Modeling 
6.1. Stellar Mass 

We use the SED fitting procedure from [Zhang] 



et al. (2012) to model the stellar mass and our re- 



sults are presented in Figure [20] Briefly speaking, 
the Hubble time is logarithmically divided into six 
independent age bins, and the star formation rate 
(SFR) in each age bin is assumed to be constant. 
A library consisting of ~ 4 x 10 6 different star 
formation histories (SFHs) is created by changing 
the relative star formation intensity among the six 
age bins. In addition, the metallicities and inter- 
nal dust extinction are also free parameters when 
creating the library. The usually adopted two- 
component (e.g. exponential plus bursts) spec- 
tral energy distribution (SED) modeling method 
can introduce a bias in the estimation of param- 
eters (Zhang et al. 2012). We use Johnson's U- 
, B-, and V-band, Ha, 2MASS J-, H-, and Ks- 
band, and Spitzer 3.6 fim photometry data in the 
SED fit. In the inner regions (radii < 1.3 kpc) 
of NGC 1569, we do not use the 3.6 /im data be- 
cause the starburst emission in these inner regions 
is expected to have significant contributions from 
hot dust (Zhang et al. 20101. In the outer re- 



properly model the outer disk. Instead of looking 
for the best-fitting parameter, we use the Bayesian 
technique to find the most probable estimates of 
parameters by constructing a probability density 
function and the corresponding cumulative distri- 
bution function for each parameter, and the con- 
fidence interval is defined as the central 68% of 
the cumulative distribution function. Any possi- 
ble degeneracies between different parameters are 
reflected in the confidence intervals. 

The right panel of Figure [20] shows the cumula- 
tive summed stellar mass as a function of radius. 
The profile approaches the horizontal asymptote 
at a mass of 2.8 x 10 8 M Q , giving the total stellar 
mass of NGC 1569. The left panel of Figure [20] 
shows the surface density mass profile, corrected 
for inclination. We see that the stellar mass den- 
sity profile follows a broken or double exponential. 
For these data, we make no assumptions about 
disk scale height or constant mass-to-light ratio. 
When we compare a de Vaucouleur's R 1 / 4 profile 
to the surface mass density in Figure [20] it is clear 
that it does not fit the data as well as a double 
exponential profile. 

To see how the R 1 / 4 profile compares to the 
more traditional surface brightness profile, we plot 
our 3.6 /im photometry surface brightness in Fig- 
ure [21] We observe a marginally better fit by the 
double exponential profile over the R 1 / 4 profile. 
We see that the R 1 / 4 profile overestimates the sur- 
face brightness in the center and underestimates it 
at ~0.5 kpc. We conclude that the surface bright- 
ness, as well as mass profile of NGC 1569 is best 
fit as a double or broken exponential disk. 

In Figure |21| there is a break in the exponen- 
tial profile around 1.3 kpc where it becomes flatter. 
This surface brightness profile is a Type III "an- 
titruncated" profile, as described by Erwin et al. 
(2005). A Type III profile is ascribed to disk flar- 
ing in the outer disk in M94 by |Herrmann fc Ciar- 



dullo (2009). They attribute this characteristic to 



gions (radii > 1.3 kpc), we use all data except 
the 2MASS data because they are too shallow to 



a recent interaction or merger. If this type of sur- 
face brightness profile is generally indicative of a 
flared disk resulting from an interaction or merger, 
then this would be additional evidence for an out- 
side influence on NGC 1569. 

To check the validity of our derived stellar sur- 
face mass density profile, we turn to our observed 
stellar kinematics. From Swaters ( 1999 and ref- 



erences therein), we use the relationship between 
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vertical stellar velocity dispersion, a z , and stellar 
surface mass density, 

g z = v/^Gz E(i?) (5) 

Here, z Q is the vertical scale-height, which we as- 
sume is constant in this relationship and is defined 
as twice the exponential scale height, z Q = 2h z . 
We also assume a vertical sech 2 scale-height distri- 
bution for the stellar component with the empiri- 



cal relationship, h/z a = 5 ( |van der Kruit fc Searle 
1981), where h and z are the disk scale- length 



and scale-height, respectively. From an exponen- 
tial disk fit to the inner part of the 3.6 /mi surface 
brightness profile, we derive a scale-length, h, of 
1.54 kpc, giving a scale-height, z Q , of 0.31 kpc for 
NGC 1569. 

Figure [22] shows our observed stellar velocity 
dispersions, plotted so that all radii are positive, 
along with the model velocity dispersion profile 
derived from Equation [5j In Figure [22] there is 
general agreement between the observed and pre- 
dicted stellar velocity dispersions in the center, 
however, the observed stellar velocity dispersions 
appear to increase as a function of radius, which 
is opposite of what is expected. This anomaly can 
be explained as an isolated area where the stars 
have a higher dispersion because they are the stars 
formed in the expanding shell described in Section 



4.2 These high dispersion stellar velocities coin- 
cide with the the peculiar dip, identified by the 
black arrow, at the edge of the expanding shell in 
Figure [TI] (R = -27" 0.5 kpc). We do not 
know how the stars act beyond ^0.5 kpc, as our 
spectra do not probe past this radius, so, we can- 
not definitively determine what happens beyond 
this increase. 

We used the stellar surface density profile, de- 
termined from the SED fitting, to model the ex- 
pected rotation velocities of the stars (rotmod), 
assuming that the total mass of the galaxy comes 
from the stars. The results for this model, cor- 
rected for inclination, are shown by the asterisk 
symbols in Figure [IT] This graph shows that in 
the inner 1 kpc, half of the observable stellar disk 
of NGC 1569, the gravitational potential is domi- 
nated by the observed stellar mass. 

6.2. Gas Mass 

The total H I mass of NGC 1569, summed over 
the integrated H I map, Figure [5] is 2.3 x 10 s M . 



To calculate the expected velocity of the gas as 
a function of radius, we use the H I mass surface 
density profile. We assume an axisymmetric, in- 
finitely thin disk comprised entirely of gas and an 
inclination of 69° ± 7° (see |Oh et al.|[20PT] for 



more details) . The model H I velocities as a func- 
tion of radius are shown by the triangle symbols in 
Figure |17| The maximum rotation velocity of the 
model H I rotation curve is 18 km s _1 , less than 
half of the asymmetric drift corrected Vm ax (filled 
circles). 

NGC 1569 also contains molecular hydrogen 
(H 2 ). According to Israel (1988), the amount of H 2 
in NGC 1569 as a whole is about 5 x 10 7 M cor- 
rected to a 3.36 Mpc distance. This produces the 
ratio Mh 2 /Mhi = 0.2 over the galaxy. The esti- 
mate of the total H2 is an approximation, though, 
because Israel determined the total H2 mass by 



scaling the results derived by Young et al. (19841 



for the central 50" of the galaxy. Also, the CO to 
H 2 conversion factor used by Israel is potentially 
dependent on metallicity, and NGC 1569 has a low 
metallicity. However, the scaled mass agrees with 



Obreschkow & Rawlings (20091 who find an H 2 
mass of ~ 2 x 10 Y M , corrected for a distance of 
3.36 Mpc. In either case, these estimates suggest 
that the total amount of H 2 is on the order of 10% 
of the total baryonic mass (stars + gas) in NGC 
1569. 

6.3. Dark Matter Mass 

In the previous sections, we have modeled the 
contributions to the rotation curve of NGC 1569 
due to stars and gas. Here, we explore the dark 
matter content in NGC 1569 by assuming that the 
observed, asymmetric drift corrected, rotation rep- 
resents motion under the total gravitational poten- 
tial of the system. Using our Vmax = 50 km s _1 , 
we calculate a dynamical mass, Afd yn , of 1.3 x 10 9 
Mq within a radius of 2.2 kpc, using the following 
equation: 

M dyn (R) = (6) 

where, V maK is the maximum rotation speed de- 
termined from the asymmetric drift corrected ro- 
tation curve, and R is the radius where V max is 
measured. We determine an H I mass and a stel- 
lar mass inside a radius of 2.2 kpc of 1.7 x 10 s M 
and 2.5 x 10 8 M , respectively. To account for 
helium, we scale the H I mass by a factor of 1.36 
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( de Blok et al.|2008 ), and achieve a total baryonic 
mass of 4.8 x 10 s M©. Subtracting this from the 
dynamical mass, we are left with only 8.2 x 10 8 
M Q of missing, dark matter within a radius of 2.2 
kpc. 

There are, however, a number of uncertainties 
in the dark matter determination. The most im- 
portant uncertainty is that we assume that we 
have found the turnover radius and, hence, the 
maximum rotation, V max , of NGC 1569. However, 
the signal-to-noise is low in the region where the 
turnover radius occurs. If the true maximum ro- 
tation speed of the galaxy is actually higher, then 
our calculation of dynamical mass, along with the 
deduced dark matter mass in the central region 
of NGC 1569, will be underestimated. Another 
uncertainty comes from the fact that these calcu- 
lations use a Vmax derived under the assumption 
that the inclination of the H I disk is 69° ± 7° , as 
determined from the tilted ring model. If a Vmax 
of 50 km s _1 is correct, then NGC 1569 has a low 
dark matter content inside the visible disk, which 
is unlike most other dwarf galaxies studied (e.g. 
Oh et al.|[20TT) |Sawala et all[20TTj [Walker et al.| 
20091). 



To see the effect of such a low dark matter 
content, we examine the dark matter fraction, 
7dm ( m^)' as a function of radius using the fol- 



lowing equation given in Oh et al. (2011 1: 



7dm (R) 



M DM V Iot (R) 2 -V star (R) 2 -V gas (R) 2 



M t , 



V TOt (Rf 



(7) 

where V rot is the rotation velocity corrected for 
asymmetric drift derived from H I, I4tar is the stel- 
lar rotation velocity assuming all the mass is in the 
stars determined from the SED fitting, and V^ as 
is the gas rotation velocity determined from the 
H I surface density assuming Newtonian dynam- 
ics. Figure [23] shows the fraction of dark matter 
function of radius. The average values in the 
inner 1 kpc are zero, within the errors, due to the 
dominance of the stars. However, in the outer re- 
gion, beyond 1.25 kpc, the fraction of dark matter 
has an average of ~0.5, which is below the average 
of 0.7 found for a sample of seven dwarf galaxies by 
Oh et al. We see here that the stars account for all 
of the dynamical mass within a radius of ~1 kpc - 
no dark matter is needed within this radius. Thus, 
NGC 1569 appears to be dark matter deficient in 



the inner disk compared to other dwarfs. 

In ACDM, the dark matter halos hosting galax- 
ies have an approximately universal shape that fol- 



lows a concentration-mass relation (e.g. Bullock et 

Hence, CDM halos 



al. 


2001 


Prada et al. 


2012) 



are a one-parameter family, typically parameter- 
ized by mass within the virialized extent (M™ ) 
or maximum circular velocity (V, nax ). Adding the 
assumption that there is a monotonic mapping of 
dark matter halo mass to galaxy luminosity (e.g., 
the brightest galaxies are hosted by the most mas- 
sive dark matter halos), a galaxy's dark matter 
halo can be uniquely determined from its lumi- 
nosity using the technique of abundance matching 



(Kravtsov et al.||2004 


Vale & Ostriker|2004 


Con- 


roy & Wechsler 


2009 


. Here, abundance match- 



ing is performed, as described in detail in Tollerud 



et al. (2011a I, using the halo V max function from 



the Millenium II simulation (Boylan-Kolchin et 
al. 2009 1 and the r band galaxy luminosity func- 
tions of Blanton et al. ( 2005 1 . The r band ab- 
solute magnitude for NGC 1569 (M r = -17.6) 
is determined from the SED fit described in Sec- 
tion |6.1| and this is used to infer that the halo of 
NGC 1569 has V max = 95 km s" 1 . This implies 
M v j r Ri 1.6 x lO n M0, and the virial extent of the 
halo is i? v i r 140 kpc. 

With these values for the total mass of NGC 
1569 's dark matter halo, the concentration-mass 
relation implies a unique shape for the halo pro- 
file. As outlined in |Tollerud etal.|(2011b), the halo 



shape is parameterized using a Navarro, Frenk, 
White| ( |1997[ ) (NFW) halo profile and the shape 
parameters are assigned so that the halo falls on 
the median of the concentration-mass relation. 
Figure [24] shows the NFW density profile as a 
function of radius for NGC 1569. The resulting 
circular velocity profile is given as the solid line 
in Figure [25] Figure [IT] shows the inner region of 
the NFW profile, solid line, as compared to the 
modeled gas and stellar masses. 

As is clear from Figure [17] the stellar mass is 
close to the expected dark matter mass in the inner 
1 kpc, and the uncertain effect of the baryons on 
the dark matter halo should be significant here. 
Thus, it is ambiguous whether or not the domi- 
nance of the baryons in the inner 1 kpc implies 
a halo profile inconsistent with ACDM, or simply 
that the baryonic mass of the galaxy is altering the 
shape of the dark matter halo in the region where 
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the dark matter does not dominate the gravita- 
tional forces. 

By R ~ 2 kpc, however, the ACDM prediction 
lies within the error bars of the measured rotation 
curve. If the model is correct, the hint of a turn 
over in the rotation curve that we see is likely to be 
an artifact of low signal-to-noise, as the rotation 
curve implied by this halo profile does not peak 
until R s=s 20 kpc, far beyond the baryonic extent 
of the galaxy. Additionally, it implies that in a 
global sense, NGC 1569 is not dark matter defi- 
cient because at a radius of ~2 kpc, the rotation 
curve is consistent with ACDM predictions. 

7. Summary 

7.1. Stellar Kinematics 

1) A systemic velocity, V^ ys , of -82 ± 7 km s^ 1 was 
determined from averaging the heliocentric radial 
velocities of the stars along the minor axis. 

2) NGC 1569 has a stellar kinematic major axis 
equal to the morphological major axis at a posi- 
tion angle, measured east from north, of 121°. 

3) The stars and gas kinematically follow each 
other. 

4) NGC 1569 has an average stellar velocity dis- 
persion of < er z > = 21 ± 4 km s _1 and a kine- 
matic measure, V ma , x /a z = 2.4 ± 0.7, which is sug- 
gestive of a thick disk. 

5) The stellar velocity dispersion increases from 21 
± 4 km s _1 to 33 ± 9 km s" 1 in the region of the 
SSCs at a radius of <~ +0.5 kpc from the galaxy 
center, which is consistent with velocities of stars 
that formed in shocked gas in the expanding shell 
around SSC A. The region in the vicinity of the 
SSCs is a disturbed part of the disk. 

6) NGC 1569 follows the observed trend of a z - 
Mb in large spiral galaxies (Swaters 1999). 

7) The 3.6 light profile of NGC 1569 follows 
a double exponential disk better than a de Vau- 
couleurs's R 1 / 4 profile. The light profile is shal- 
lower in the outer disk than in the inner disk. 

7.2. Ionized Gas Kinematics 

1) Multiple velocity components are observed in 
the spectra as a result of mechanical energy input 
from the star burst. 

2) There is one expanding shell observed from the 
high velocity ionized gas seen along the major axis 
and indications for outflow. 



3) The stellar velocities match low velocity ionized 
and H I gas and are likely produced by the stars 
that caused the expanding shell (high velocity gas) 
observed on the western side of the disk. 

7.3. H I Kinematics 

1) Extended H I emission to the south and north 
of NGC 1569 was imaged for the first time. 

2) NGC 1569 shows a dominance of non-circular 
motions across the H I disk as seen in the intensity- 
weighted velocity field. 

3) A non-circular motion (NCM) H I cloud was ob- 
served and appears to be impacting NGC 1569 at 
the eastern edge of a dense H I cloud that lies west 
of the SSCs, as determined from a double Gaus- 
sian decomposition method that was successfully 
applied to our H I data. 

4) We suggest that the impact from the NCM H I 
cloud caused the dense ridge in the center of the 
disk to become unstable, thus creating the star- 
burst and enabling the formation of the SSCs. 

7.4. Kinematic and Morphological Com- 
parisons 

1) The disturbed stellar kinematics correspond to 
the region west of SSC A. 

2) A "butterfly" shape in the H I gas morphol- 
ogy is observed, likely caused from outflow from 
the starburst. Gas is escaping from the body of 
the butterfly, perpendicular to the stellar and gas 
disks. 

7.5. Mass Modeling 

1) A total stellar mass of 2.8 x 10 8 M Q was deter- 
mined from a new SED fit and a total gas mass, 
including helium and metals, of 2.3 x 10 8 M was 
determined. 

2) A dynamical mass of 1.1 x 10 9 M Q was calcu- 
lated from the V maK = 50 km s" 1 at a radius of 
2.2 kpc. 

3) The stars dominate the gravitational potential 
in the inner 1 kpc of the H I disk. 

4) The derived rotation curve from the H I kine- 
matics compared to the modeled H I rotation from 
the surface mass density and the modeled stellar 
rotation, shows that only a small amount, 80% of 
the baryonic mass, is needed in dark matter. This 
is unlike typical dwarf galaxies. 

5) The rotation of the stars and gas seem to agree 
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with an NFW profile at a radius of <~ 2 kpc. 
6) The NFW model predicts a virial mass of 1.6 x 
10 11 Mq and gives a virial radius of 22 kpc, well 
beyond the radius where baryonic material is de- 
tected. 
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Fig. 1. — Outline of all major gas features in NGC 1569. The black contours in the left image are from the 
integrated H I data and the center of the map is left blank so that the scale of the HST image and NCM 
H I cloud and UD H I cloud can be seen. The extended H I to the northeast and south of NGC 1569 were 
imaged for the first time in this work. The right image is an HST composite image showing the Ha features 
along with the UD and NCM H I clouds. SSCs A and B are identified in the right HST image as well as the 
position, location and size of the optical slit used to derive the stellar kinematics. The white curved arrows 
in both images show the rotation of the disk. The near and far sides of the disk have been labeled as well 
as the relative velocities of the H I clouds and features. Credit: NASA, ESA, the Hubble Heritage Team 
(STScI/AURA), and A. Aloisi (STScI/ESA) 
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Fig. 2. — HST image from Figure [T] showing the slit positions and the location of SSCs A and B. The bar 
at the bottom left in both panels shows ~ 1 kpc in distance and the blue star marks the kinematic center 
determined from the H I tilted ring model. 
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Fig. 3. — Top: Five 1800s spectra co-added to create a two-dimensional galaxy spectrum of the major axis 
of NGC 1569 before sky subtraction. Bottom : Same as top only after sky subtraction. For both panels, 
wavelength is along the horizontal axis (~238 A total coverage) and the vertical axis is the cross-dispersion 
direction along the slit (~3' ). The right panel shows residual sky emission lines, indicated by the arrows. 
Only the bright region in the spatial center was used in the cross-correlation process. 
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Fig. 4. — Top: Co-added, sky subtracted, two-dimensional spectrum of the O III A5007 A emission feature 
for the major axis of NGC 1569. Bottom : Same as left, only for the minor axis. Wavelength is along the 
horizontal axis (~41 A total coverage) and the vertical axis is the cross-dispersion direction of the slit (~3' ). 
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Fig. 5. — Integrated H I flux of NGC 1569. The multi-scale cleaning algorithm implemented by LITTLE 
THINGS was able to image tenuous emission to the south and northeast for the first time. 
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Fig. 6. — Same HST image as in Figure [T] only shown in greyscale with the major axis slit superposed. 
The ovals are the size and position of the extracted stellar spectra. The red ovals correspond to red-shifted 
velocities relative to the systemic velocity whereas blue ovals correspond to blue-shifted velocities. The green 
'X' marks the center of the slit and the galaxy. SSC A and B are marked by the green circles. Credit: NASA, 
ESA, the Hubble Heritage Team (STScI/AURA), and A. Aloisi (STScI/ESA) 
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Fig. 7. — Example of three, one-dimensional galaxy spectra (first three panels) and one template spectrum (last 
panel) for comparison. The galaxy spectra, observed along the major axis of NGC 1569, are representative of the 
range of signal-to-noise (S/N) and span the spatial area of the slit; the top panel has a S/N of ~ 7 at a radius of 
+18" (east) from the center, the second panel has S/N of ~ 24 at a radius of -4" (west) from the center, the third 
panel has a S/N of ~ 6 at a radius of -34" (west) from center. The stellar template spectra has a S/N of ~ 100, 
which was typical throughout both observing runs. We used the regions shown by the horizontal brackets in the 
cross-correlation method (CCM). The Mg lb absorption features are marked in each panel. 
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Fig. 8. — A single extracted spectrum of the O III A5007 A emission feature at a radius of -13" from the 
center of NGC 1569 along the major axis. Three bright Gaussian peaks are seen. The total width of the 
blended emission feature at zero intensity is ^4.1 A. or, ^176 km s _1 . 
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Fig. 9. — Intensity-weighted H I velocity field of NGC 1569. The red outline identifies the region over which 
the bulk motion extraction procedure was performed. The black line marks a ~1 kpc distance scale. The 
tiny dot in the lower right corner near (04:29:35, +64:35:00) shows the resolution. 
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Fig. 10. — Results of the kinematic deconvolution procedure of Oh et al. (2008): (a) Bulk rotation map. 
(b) Strong non-circular motion map. (c) Weak non-circular motion map. (d) Combined strong + weak 
motion map. The black oval identifies a strong non-circular motion (NCM) H I cloud at -50 km s _1 , clearly 
redshifted with respect to the underlying (bulk) motion. The black line in the center of each velocity field 
identifies the position and length of the slit along the major axis used in the stellar spectroscopy. The black 
line in the center of the figure shows approximately 1 kpc distance. 
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Fig. 11. — Observed heliocentric velocities in the stars along four PAs. Red diamonds, green stars, and 
black filled circles represent independent points derived from different observing nights. The dashed lines 
represent the ±6.9 km s" 1 uncertainty in V sys . East are positive distances from the galaxy center. Error 
bars along the x-axis show the angular area summed along the slit. 
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Fig. 12. — Position-velocity diagrams taken along the four optical slit PAs showing H I contours with stellar 
data over plotted. Positive distance corresponds with eastward direction. The solid line in the upper left plot, 
major axis, shows the ordered rotation of the H I. The oval in the upper left plot identifies the "forbidden" 
H I motion. The stars and the H I follow each other closely at all four PAs. 
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Fig. 13. — Stellar velocity dispersions observed along the four optical slit PAs. Black filled circles, red 
diamonds, and green stars represent independent points derived from different observing nights. The solid 
line shows the average velocity dispersion, < er z > = 21, and the dotted lines indicate the uncertainty, ± 
4.2 km s _1 . Error bars along the x-axis show the angular area summed along the slit. East are positive 
distances from the galaxy center. 
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Fig. 14. — Position- velocity plot of the ionized gas (black filled circles), H I gas (blue boxes), and stars (yellow 
stars) along the major axis. The solid red outlines an expanding gas shell. Note the close correspondence 
between the stars and ionized gas, including peculiar motions, in the western-most oval as indicated by the 
black arrow. East are positive distances from the galaxy center. 
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Fig. 15. — Measured rotation curve of the H I gas determined from the bulk motion map using a tilted 
ring model. A V ro t of 33 km s _1 is determined from averaging over the region where the curve is fiat from 
~2.0 - 3.2 kpc. The errors on the rotation curve points are the standard deviations calculated by assuming 
a solid body rotation out to 2.2 kpc. The average systematic error, determined from all observational and 
computational uncertainty, is shown for reference. 
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Fig. 16. — Velocity dispersion map of the bulk velocity component only for all extracted velocity peaks above 
2a. 
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Fig. 17. — Mass modeling of both the stars (asterisks) and the gas (triangles) along with the observed bulk 
rotation curve (boxes) and the asymmetric drift corrected rotation curve (filled circles). The NFW profile 
(solid line) is also shown. The average uncertainty for the model stellar rotation are shown in the lower right 
corner for comparison to the asymmetric drift corrected rotation curve. Within the uncertainties, the stellar 
mass dominates the entire rotation curve out to 1 kpc, half of the observable stellar and gas disks. 
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Fig. 18. — We show our deep Ha image with the integrated H I contours over plotted. More Ha streamers 
and arcs can be seen to the south due to the lower H I density. To the north, we observe a strong Ha 
chimney that extends all the way through the H I disk. The ultra-dense (UD) H I cloud is identified as well 
as the two dense clouds on the opposite side of the optical galaxy. We believe these features were at one 
time connected and extended through the current optical galaxy. The non-circular motion (NCM) cloud, 
outlined by the blue oval, may have caused the dense H I gas to collapse into the starburst that ended nearly 
20 million years ago. 
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Fig. 19. — Intensity-weighted velocity contour map over-plotted on our Ha image. The features identified by 
Stil & Israel (1998) are marked, as well as the position of the NCM cloud shown in Figure 10 We show two 
new features here that make up what we call butterfly wings. The cavity between the two wings matches 
the position of X-ray emission ( Martin et al.|2002 l. We also show the position of the UD H I cloud. The Ha 
tail is almost completely within the NCM H I cloud but has the same velocities as the UD cloud. Therefore, 
we conclude that this tail is the lit up edge of the UD H I cloud. 
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Fig. 20. — Left : Logarithmic stellar mass surface density profile from fitting the SED in elliptical annuli. 
This has been corrected for an inclination of 69° ± 7° . The solid lines are a fit to the profile and represent a 
double exponential disk. The dashed line represents an R 1 / 4 fit to the surface density profile where we have 
assumed a constant mass-to- light ratio and constant disk scale height. An exponential disk is a better fit 
than an R 1 / 4 . Right: The logarithmic integrated stellar mass for NGC 1569 from the elliptical annuli SED 
fitting. 
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Fig. 21. — Inclination corrected 3.6 /zm surface brightness profile. The solid lines are an exponential fit to 
the profile, and the dashed line represents an R 1 / 4 fit. The broken exponential profile fits the data better 
than the R 1//4 profile. We show this figure because it is the traditional way that surface brightness profiles 
are studied and the transformation to mass surface density, shown in Figure |20| is not linear. 



41 




Fig. 22. — Our observed stellar velocity dispersions measured from our spectra (star symbols). The solid line 
represents the inclination-corrected model velocity dispersion expected from the stellar mass surface density, 
determined from SED fitting. z is the stellar scale height. The bulk of the points (R < 0.35 kpc) lie on or 
near the expected values. But, the points between 0.35 < R < 0.6 kpc are too high, indicating a region of 
unusually high stellar velocity dispersion. We interpret this as a disturbed region where the stars formed in 
an expanding shell. 
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Fig. 23. — Dark matter fraction, "/dm, as a function of radius. The negative dark matter fractions in the 
inner 1 kpc are due to the dominance of the stellar mass to the rotation curve. Within the errors, the stars 
account for all of the mass within a radius of 1 kpc. Other dwarfs, by contrast, have -j dm ~ 0.7 on average 
over their entire disks. 
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Fig. 24. — NFW density profile from the abundance matching method, assuming an M r — —17.6. 
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Fig. 25. — NFW velocity profile for the abundance matching method, assuming an M r = —17.6. 
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Table 1 

Global Parameters for NGC 1569 



Parameter 


Value 


Ref 


Other Names 


UGC 03056, ARP 210, 


1 




VII Zw 016 




D (Mpc) 


3.36 ± 0.20 


2 


M v 


-18.2 


3 


H25 (mag arcsec~ 2 ) 


22.3 


5 


Galaxy diameter to 25 mag arcsec -2 in B, D25 (arcmin) 


3.6 


3 


V-band disk scale length, Kd (arcmin) 


0.39 ± 0.02 


3 


Center (RA, DEC) (J2000) 


(04:30:49.8, +64:50:51) 


3 


Minor-to-major axis ratio, b/a 


0.55 


3 


Lira (ergs s" 1 ) 


5.7 x 10 40 


4 


Star Formation Rate (SFRo) (Moyr^kpc -2 ) 


1.3 


4 


H 2 Mass (M e ) 


5 x 10 7 


6 


Stellar parameters determined from this work 


V sys (km s -1 ) 


-82 ± 7 




Major Axis PA (degrees) 


121 




Inclination from de-projection, i op t (degrees) 


60 ± 10 




Average Velocity Dispersion, < cr z > (km s _1 ) 


21 ± 4 




Velocity Dispersion near SSC A, <r z (km s -1 ) 


33 ± 9 




Total Stellar Mass (M©) 


2.8 x 10 8 




Optical scale-length, h (kpc) 


1.54 




Optical scale-height, z (kpc) 


0.31 




H I parameters determined from this work 


V sys (km s 1 ), 


-85 




Major Axis PA (degrees) 


122 




Inclination from tilted ring model, iui (degrees) 


69 ± 7 





Kinematic Center (X pos , Y poB ) (RA, DEC) 
Kot (km s _1 ) 

Vmax (km S _1 ) 

Average Velocity Dispersion, < <thi > (km s _1 ) 
Total Gas Mass (M ) 

Mdyn (M Q ) 

Mom (M ) 
M vir (M ) 
7? vir (kpc) 

Knax z 



(04:30:46.125, +64:51:10.25) 
33 ± 10 
50 ± 10 
10 

2.3 x 10 s 

1.3 x 10 9 
8.2 x 10 8 
1.6 x lOH 

140 

2.4 ± 0.7 



References. — (1) NASA Extragalactic Database; (2) Grocholski et al. ( 2008 1;(3) Hunter fc 



Elmegreen (2006 1; (4) Hunter fc Elmegreen (2004 1; (5) de Vaucouleurs et al. (1991 1;(6) Israel (1988) 
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Table 2 

Optical spectroscopic observations of NGC 1569 taken with the KPNO 4-m + Echelle 

SPECTROGRAPH IN 2008. 



UTDate PA a iV gal b T gal (s) c N sky d T sky (s)° iV spcc f AW 

Obs. 



Mg lb spectral region 



02/01 


Major axis (121°) 


5 


9000 


5 


9000 


7 


4 


02/02 


Major axis (121°) 


5 


9000 


5 


7393 


9 


5 + 1 


02/03 


Minor axis (31°) 


4 


7200 


4 


5673 


4 


11 + 1 


11/19 


+ 45° from major axis (166°) 


8 


13,800 


7 


12,000 


4 


5 


11/20 


— 45° from major axis (76°) 


3 


5400 


2 


3600 


2 


6 + 1 


11/21 


— 45° from major axis (76°) 


8 


14,400 


8 


13,860 


6 


9 


11/22 


— 45° from major axis (76°) 


5 


9000 


5 


8841 


3 


10 




Major axis (121°) 


3 


5400 


3 


5400 


4 


10 


O III A5007 A observations 


11/23 


Major axis (121°) 


3 


4800 


3 


1800 


95 


NA 


11/23 


Minor axis (31°) 


3 


1800 


1 


600 




NA 



a Position angle observed; 

b Number of galaxy exposures co-added into a single spectrum; 

Total integration time, in seconds, on galaxy; 

d Number of sky exposures co-added into a single spectrum; 

"Total integration time, in seconds, on sky; 

f Number of one-dimensional spectra extracted; 

g Number of radial velocity standard stars used in cross-correlation function. The is in 

reference to the solar spectrum, from twilight flats, used as a radial velocity standard star. 
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Table 3 

Radial velocity standard stars used for cross-correlation. 



HD 


Apparent 


Metallicity 


Spectral 


RA (J2000) 


DEC (J2000) 


Vhelio a CTV 


T exp UT Date 


No. 


Magnitude, V 


Fe H 


Type & Class 


hh:mm:ss 


. / // 


km s _1 km s _1 


(sec) Obs. (2008) 



8779 


6.41 


-0.40 2 


K0 IV 


01:26:53.5 


-00:21:18 


-5.0 


0.5 


36 


02/03 


9138 


4.84 


-0.37 3 


K4 III 


01:30:37.9 


+06:11:15 


+35.4 


0.5 


7 


02/03 


















10 


11/22 


12029 


7.44 




K2 III 


01:59:11.2 


+29:25:16 


+38.6 


0.5 


90 


11/22 


22484 


4.28 


-0.10 3 


F9 IV-V 


03:37:18.5 


+00:25:41 


+27.9 


0.1 


5 


02/01 


















5 


02/02 


















5 


02/03 


















5 


11/22 


23169 


8.50 




G2 V 


03:44:23.9 


+25:45:06 


+13.3 


0.2 


240 


02/03 


26162 


5.50 


-0.02 1 


Kl III 


04:09:39.8 


+ 19:37:52 


+23.9 


0.6 


15 


02/01 


32963 


7.60 


0.08 1 


G5 IV 


05:08:27.3 


+26:20:18 


-63.1 


0.4 


80 


02/02 


















80 


11/21 


65583 


6.97 


-0.56 3 


G8 V 


08:01:03.6 


+29:11:09 


+12.5 


0.4 


60 


11/19 


















75 


11/22 


65934 


7.70 




G8 III 


08:02:42.1 


+26:36:49 


+35.0 


0.3 


90 


11/22 


66141 


4.39 


-0.30 3 


K2 Illb Fe-0.5 


08:02:42.5 


+02:18:39 


+70.9 


0.3 


10 


11/22 


75935 


8.46 




G8 V 


08:54:20.3 


+26:52:50 


-18.9 


0.3 


170 


11/21 


90861 


6.88 




K2 III 


10:30:22.2 


+28:32:15 


+36.3 


0.4 


60 


11/21 


92588 


6.26 


-0.10 2 


Kl IV 


10:41:50.1 


-01:47:11 


+42.8 


0.1 


35 


11/19 


















40 


11/22 


102494 


7.48 


-0.26 4 


G9 IVw... 


11:48:22.8 


+27:17:36 


-22.9 


0.3 


90 


11/22 


122693 


8.11 




F8 V 


14:03:15.6 


+24:31:14 


-6.3 


0.2 


170 


02/03 


126053 


6.27 


-0.45 3 


Gl V 


14:23:41.4 


+01:12:08 


-18.5 


0.4 


35 


02/02 


















35 


02/03 


132737 


7.64 




K0 III 


15:00:14.4 


+27:07:37 


-24.1 


0.3 


80 


02/01 


















80 


02/02 


















80 


02/03 


136202 


5.06 


-0.08 3 


F8 III-IV 


15:19:44.9 


+01:44:01 


+53.5 


0.2 


24 


02/03 


144579 


6.66 


-0.69 1 


G8 IV 


16:05:14.4 


+39:08:02 


-60.0 


0.3 


35 


02/01 


















35 


02/03 


145001 


5.00 


-0.26 5 


G5 III 


16:08:27.6 


+ 17:01:29 


-9.5 


0.2 


10 


02/02 






-0.04 6 












10 


02/03 


154417 


6.01 


F8.5 IV-V 


17:05:42.8 


+00:41:26 


-17.4 


0.3 


20 


02/03 


182572 


5.16 


-0.37 1 


G7 IV R8 1 


19:25:22.5 


+ 11:57:47 


-100.5 


0.4 


24 


11/20 


















24 


11/21 


187691 


5.11 


0.07 7 


F8 V 


19:51:26.1 


+ 10:26:15 


+0.1 


0.3 


24 


11/19 


















24 


11/21 


194071 


7.80 




G8 III 


20:22:58.8 


+28:16:26 


-9.8 


0.1 


90 


11/20 


203638 


5.41 


0.02 5 


K0 III 


21:24:38.4 


-20:48:55 


+21.9 


0.1 


15 


11/20 


















24 


11/21 


212493 


4.79 




K0 III 


22:28:17.3 


+04:44:19 


+54.3 


0.3 


7 


11/19 


















7 


11/21 


















10 


11/22 


213014 


7.45 




G9 III 


22:28:36.3 


+ 17:18:25 


-39.7 


0.0 


75 


11/19 
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Table 3 — Continued 



HD 


Apparent 


Metallicity 


Spectral 


RA (J2000) 


DEC (J2000) 


Viielio a 


(TV 


T cx p 


UT Date 


No. 


Magnitude, V 


Fe H 


Type & Class 


hh:mm:ss 


O . / , // 


km s _1 


km s" 1 


(sec) 


Obs. (2008) 


213947 


6.88 




K2 


22:35:00.6 


+26:38:32 


+ 16.7 


0.3 


40 


11/21 


222368 


4.13 


-0.15 3 


F7 V 


23:40:23.3 


+05:40:21 


+5.3 


0.2 


5 


11/20 


















5 


11/22 


223311 


6.07 




K4 III 


23:48:58.7 


-06:20:00 


-20.4 


0.1 


20 


11/20 



30 11/21 



a Heliocentric radial velocities. 

All data in table, except for metallicity and the last column, come from the United States Nautical Almanac (2008). Metallicities 
are listed, where available, and references are given below. 



1 Soubiran et al. 



Randich et al. 



( 2008 1 
( 1999 1 
(20071 



lYoss fe Griffin] ( |1997 1 



' Cenarro et al 



McWilliam 



£ Chen et al. 



Fuhrmann 



(19901 
( 2000 1 
(19981 
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